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Abstract 

We introduce a Multifractal Random Walk (MRW) defined as a stochastic integral 
of an infinitely divisible noise with respect to a dependent fractional Brownian motion. 
Using the techniques of the Malliavin calculus, we study the existence of this object 
and its properties. We then propose a continuous time model in finance that captures 
the main properties observed in the empirical data, including the leverage effect. We 
illustrate our result by numerical simulations. 
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1 Introduction 

Starting with the seminal work of Mandelbrot about cotton price [21], several studies of 
financial stock prices times series, have allowed to exhibit some particularities of their 
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fluctuations. Without making a comprehensive list, we can mention the appearance in 
empirical data of the following properties: non Gaussian distributions due to now well 
known fat tails of financial returns, the so-called volatility clustering that means that the 
volatility fluctuations are of intermittent and correlated nature, scaling invariance, long run 
correlation in volatility, leverage effect and so on (see e.g [18], [14], [25] for an extensive 
review). Thus, constructing theoretical models for financial returns that include all the 
properties listed before appears as a very interesting challenge. Many scientific works, in 
economics or mathematics, proposed various models for asset returns. The ARCH model 
introduced by Engle [15] offers an interesting base of work, and after the seminar work by 
Engle a vast literature on ARCH and related models has been developed. One of the first 
extensions of the ARCH model, called the GARCH model, has been introduced by Bollerslev 
in [9] and it also has been the object of various generalizations. While the GARCH models 
are capable of capturing volatility clustering, there are some drawbacks of the model. For 
example, the GARCH models are unable to represent volatility asymmetry. Due to the 
presence of the squared observed data in the conditional variance equation, the positive and 
negative values of the lagged innovations have the same effect on the conditional variance. 
In the finance literature, it has been noticed that volatility often is affected by negative 
and positive shocks in different ways. Another inconvenient of the GARCH model is the 
fact that, to ensure positiveness of the conditional variance, non-negative constraints on the 
coefficients in the variance equation must be imposed. To take into account the asymmetric 
effects on conditional second moments and to avoid the non-negativity constraints on the 
coefficients in the variance equation, in [23] the so-called EGARCH model has been proposed 
by D. Nelson. 

On the other hand, many recent empirical studies, based on huge data sets available 
nowadays, put in light new aspects of the financial returns. For example, they suggested 
that the fluctuations of the asset process displays multifractal properties (see e.g. [16] or 
[17]). Taking into account these multifractal character, several authors proposed models 
based on the on "cascade" random processes and "multifractal random walk (MRW)". We 
refer, among others, to [12], [5] or [4]. Usually, the noise in these models are defined by 



where X is a multifractal random measure and Y is a self-similar process with stationary 
increments, independent of X. Another construction defines the noise via a stochastic 
integral as (see e.g. [7], see [1]) 



where Q is a suitable fractal noise and Y is a self-similar process with stationary increments, 
independent of Q. A natural choice for the process Y, both in (1) and (2), is the fractional 
Brownian motion (fBm). 

While these multifractal models (1), (2) are able to capture several properties ob- 
served on the empirical data (scaling, volatility clustering or long-range dependance), they 



Z(t) = Y(X([0,t]), i £ M 



(1) 




(2) 
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reproduce poorly the leverage effect. The leverage effect is understood as the correlation be- 
tween the log-return at a fixed time t and his volatility (that may be defined as the squared 
or absolute log-return) varying around t, for example t — lOOAt < t' < i + lOOAi. Empirical 
studies (see e.g. [8] or [11]) have shown that this quantity is close to for past volatility and 
follow an exponential law for future volatility. This stylized fact is interpreted as a panic 
effect, the two quantities are negatively correlated just after t and go back quickly to as 
t' increases. 

We introduce here a generalization of the model (2), able to take into account the 
leverage effect and allows a more flexibility for long-range dependence in log-return. That is, 
in the model (2) we will consider the processes Q to be multifractal process called "infinitely 
divisible cascading noise" and the process Y to be a fractional Brownian motion with Hurst 
parameter H > i. The same situation has been treated in [1], [20], but in addition, we 
construct a model in which the processes Q and B H are not independent anymore. An 
alternative construction has been proposed in [2] but in this reference the integral in (2) is 
not a true stochastic integral with respect to the fBm. We will define the integral (2) using 
the techniques of the Malliavin calculus and we will study the properties of the fractional 
Multifractal Random Walk. 

Our paper is organized as follows. Section 2 contains some preliminaries on mul- 
tifractal process and the basic elements of the Malliavin calculus that we will need in our 
paper. In Section 3 we analyze the stochastic integrals that are used to define the MRW 
while in Section 4 we study the existence of the fractional MRW as a limit of a family 
of stochastic integrals. The properties of the fractional MRW (scaling, moments etc) are 
discussed in Section 5. The last section (Section 6) contains a numerical analysis of the 
data for several financial indices and we compare the simulation of our theoretical model 
with the real data. 

2 Preliminaries 

We present here he basic facts related to the infinitely divisible cascading noises and we 
introduce the basic tools of the Malliavin calculus. 

2.1 Infinitely divisible cascading noise 

Let M denote an infinitely divisible, independently scattered random measure on the set 
R x R + with generating infinitely divisible distribution G satisfying 



for some function p and for every q € M. We assume that M has control measure m on 
M x meaning that for every Borel set icRx M+ it holds 




Ee qM(A) = e -p(q)m(A) for eyery g £ R _ 
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The fact that M is independently scattered means that the random variables 

M(A 1 ),M(A 2 ),...,M(A n ) 

are independent whenever the Borel sets A±, ...,^4n £ M x M + are disjoint. We define the 
Infinitely Divisible Cascading noise (IRC) by 

e M{C r {t)) 

Qr{t) = Ee M ( c ^)) (3) 
for every r > and f G 1. Here C r (t) is the cone in K x M + defined by 

C r (t) = {(t',r'),r< r ' <l,t- r -<t' <t+ r -}. (4) 

We will use the following facts throughout our paper. We refer to [13] or [1] for the 
their proofs. First, let us note the scaling property of the moments of the IRC 

EQ r (t) q = e "P(Q)mr(0) 
and the expression of its covariance: for every r > and f,s£M 

EQ r (t)Q r (s) = e -(p(2)«v(l*--l) (5) 
where we denoted by for u > 0, r > 0. 

m r (u) = m (C r (0) n C r {u)) 

and by 

<p(q) = p(q) - ffp(i)- (6) 

The scaling of the moment of Q can be extended to the following scaling property 
in distribution : for t G (0, 1) 

(Qrt(tu)) um =W e Qt (Qr(u)) um (7) 

where " =^ " means equivalence of finite dimensional distributions. Here $7( denotes a 
random variable independent by Q, which satisfies, if the measure m is given by (14), 

Ee qQt = t qv( - q) . (8) 

Remark 1 As noticed in [1], we have </ 7 (2) < 0. 

In [1] (see also [20] ) the Multifractal random walk (MRW) based on fractional 
Brownian motion is defined as limit when r — > (in some sense) of the family of stochastic 
integrals (Z,j? (t)) r>0 defined by 

Z?(t) = [ Q r (u)dB H {u), te[0,T) (9) 
J o 
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where (B^ )te[o,T] is a fractional Brownian motion with Hurst parameter H G (0, 1). The 
fractional Brownian motion (B^ )te\o,T] with Hurst parameter H G (0, 1) is a centered 
Gaussian process starting from zero with covariance function 

In [1] it is assumed that M and -B^ are independent. Therefore the stochastic integral 
with respect to B H in (9) behaves mainly as a Wiener integral since, because of the inde- 
pendence, the integrand Q r (u) can be viewed as deterministic function for the integrator 
B H . 

Another important fact in the development of this theory is that the IRC Q is a 
martingale with respect to the argument r. Let us recall the following result (see [13]): 

Lemma 1 For every u > the stochastic process (Q r (u)) r>0 is a martingale with respect 
to its own filtration. As a consequence, for every u, v, r, r' > with r < r' it holds 

EQ r (u)Q r/ (v)=EQ r (u)Q r (v). (10) 

The property (10) plays an important role in the construction of the MRW process in [1] 
or [13]. 



2.2 Malliavin calculus 

Let (Wt)teT be a classical Wiener process on a standard Wiener space (Q, F, P). By W(tp) 
we denote the Wiener integral of the function cp £ L 2 (T) with respect to the Brownian mo- 
tion W. We denote by D the Malliavin derivative operator that acts on smooth functionals 
of the form F = g(W(<pi), . . . , W((p n )) (here g is a smooth function with compact support 
and ipi G L 2 (T) for i = l,..,n) 

n Pi 

DF = J2g^.(W( Vl ),...,W(c Pn )) ( p i . 
i=i 

The operator D can be extended to the closure D p ' 2 of smooth functionals with respect to 
the norm 



\Ff p>2 =MF 2 + 5>||D i F||| 3Ti 



i=i 



where the i th Malliavin derivative is defined iteratively. The adjoint of D is denoted 
by 5 and is called the divergence (or Skorohod) integral. Its domain {Dom{5) ) coincides 
with the class of stochastic processes u G L 2 (£l x T) such that 

\E(DF,u)\ < c\\F\\ 2 
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for all F £ D 1 ' 2 and S(u) is the element of L 2 (Q) characterized by the duality relationship 

E(F5(u)) = E{DF,u). 

For adapted integrands, the divergence integral coincides to the classical Ito integral. A 
subset of Dom{5) is the space L 1,p of the stochastic processes such that ut is Malliavin 
differentiate for every t and 



i hp n hp 
11 K „ := 11" 



+ \\Du\\ p r „ /rr2 ^ n .s < 00. 



l,p - — \\ u \\LP(Txn) ^ ll-^"!! LP (T 2 xfi) 

We will need Meyer's inequality that allows to estimate the L p moment of the Skorohod 
integral 

E\5(u)\ p <\\u\\l p . (11) 
For the L 2 moment of the Skorohod integral we have the explicit formula 

E5(u) 2 = / u 2 s ds + D r u s D s u r drds (12) 

Jt Jt Jt 

if u G L 1 ' 2 . We also recall that the Malliavin derivative satisfies the chain rule 

Df(F) = f'(F)DF (13) 
if / is a differentiable function and F £ D 1 ' 2 . 

3 The construction of the fractional Multifractal Random 
Walk with dependent noise 

Our purpose is to give a meaning to the stochastic integral (9) in the situation when the IRC 
Q and the fBm B H are not independent. We will use techniques related to the Malliavin 
calculus. In order to apply these type of techniques we will restrict ourselves to the case 
when the measure M introduced in Section 2 is Gaussian. 

3.1 The Gaussian isonormal noise 

We introduce a MRW without independence between the measure M (denote by W in 
our settings) and the integrator B H . We will restrict to the case where M is a Gaussian 
measure. More precisely, we will consider (W(h),h G H) an isonormal process, that is, a 
centered Gaussian family with 

EW(h)W(g) = (h,g) H 
for every g, h S H. The Hilbert space H will be 

H = L 2 (M x R+,B(R x M+),m)) 
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where m is the control measure. In this work we will consider 

cdv 

m(dt, dr) = dt— if < r < 1 (14) 

and m vanishes if r > 1. Here c is a strictly positive constant. This is called in [3] (see also 
[1]) the exact invariant scaling case. 

The following properties of the noise W are immediate: denote VF(^4) = W(1a) for 
A G B(R x M + ). 

• For every AeB(lxl + ) it holds that 

W{A) ~ N(0,m(A)). 

•We have 

EW(A)W{B) = 

if the Borel sets A,B C M x M+ are disjoint. This implies that the random vari- 
ables W(A) and W(B) are independent, so the random measure W is independently 
scattered. 

• For every q G R and ^4 G 0(M x E + ) we have 

which means that 

p(o) = -\<1 2 and = + I = _ 2 9 ^ ~ ^ 

Let us use the following notation: 

mi(dt) = dt and m,2(dr) = c~. 

Above mi is the Lebesque measure on R and ni2 is a measure on IR + . Clearly m = mi8m2, 
the product measure. 

For every t > and G £>(M+) such that 

m 2 (A 2 ) = 1 

(take for example ^2 = (c, 00)) we set 

W^(t):=W(l mxA2 ). (15) 
The following result is immediate. 

Proposition 1 If A2 G is such that 7712(^-2) = 1 then the process (W^\t))t>o given 

by (15) is a standard Brownian motion on the same probability space as W. 
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Proof: It is clear that is a Gaussian process. Let us compute its covariance. For 

every s, t > it holds that 

EwW(t)wM(s) = EW{l MxA2 )W(l MxA2 ) 

= ( 1 [0,t]xA 2 > 1 [0,s]xA 2 )h = (l[0,t]) 1 [0,s])L 2 (R) m 2(-42) 
= i A S 

and this implies that is a Wiener process with respect to its own filtration. ■ 

3.2 The approximating Multifractal Random Walk with dependent frac- 
tional Brownian motion 

We introduce the Multifractal Random Walk with based on the fractional Brownian motion 
by the formula 

Z?(t) = [ Q r (u)dB H (u) (16) 
Jo 

where for every r > 0, u > the integrands Q r {u) is defined by (3) with W instead of M 
and the fractional Brownian motion B is given by 



B H (t) = f K H (t,s)dW (1 \s) 
Jo 



(17) 



where K H is the usual kernel of the fractional Brownian motion and is the Brownian 
motion defined by (15). By Proposition 1, it is clear that B is a fractional Brownian 
motion. Recall that, when H > | the kernel K H (t, s) has the expression 

K H (t, s) = chs 1 / 2 - 11 f\u - s) H ~ 3 / 2 u H - l / 2 d U 

J s 

( H(2H—1) \ ^~' 2 

(see [24]) where t > s and ch = I g(o-2H h-i/2) ) anc ^ ^v> ') * s Beta function. For 



^(2-2H,H-l/2) 

t > s, the kernel's derivative is 



Ot 



it,s)=c H [-) (t-s) 



H-3/2 



In the sequel we will simply denote K H := K. 

The stochastic integral in (16) is a divergence (Skorohod) integral with respect to 
B H as defined in e.g. [24] . Actually, when H > \ we can write 

Z?(t) = J* Q r (s)dB H (s) = j* dWM{s) (J*dad 1 K(a,s)Q r (a) \ . 
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The last equality is due to the definition of the stochastic integral with respect to B (see 
[24], Chapter 5 for example). We can also express Z r (t) as a Skorohod integral with respect 
to the measure W by the formula 

Zr{t)= I [ dW(s,r')( f dad 1 K(a,s)Q r (a)) (18) 
Jo J A 2 \Js / 

where A2 is the Borel set satisfying 7712(^2) = 1 an d that appears in the definition of W^ 1 ). 

Remark 2 Clearly the fBm B and the noise W are dependent. Actually, it can be seen 
that B H and W are correlated in general. Indeed, for every Borel set iclx M + and for 
every t > 

EB H {t)W(A) = E [ K(t,s)dW i - 1) {t)W{A) 
Jo 

= E f f K(t,s)dW(s,r)W(A) 
Jo Jo 

= E I I K{t,s)ds% 

J J{(0,t)x(l,oo))nA r 

and this is in general not zero. 

Let us further discuss the dependence between the integrator B and the integrand 
Q in (9). We need to distinguish two situations. It depends on the relations of the set A2 
and the unit interval (0, 1). 

The disjoint case: We have A2 n (0, 1) = (this happens when A2 = (c, 00) with 
c > 1 for example). 

On the other hand, Bf is independent by W(C r (t)) for every s,t,r. Indeed, since 
C r (t) C R x (0, 1) (see (4)) we have 

EB?W(C r (t)) = 

and since (B H (s),W(C r (t))) is a Gaussian vector, we obtain the independence. 

The non-disjoint case: A2 n (0,1) 7^ (this happens when A2 = (c, 00) with 
c < 1). In this case EB H (s)W{C r {t)) ^ in general and so B H {s) and W(C r (t)) are 
dependent. 

We will refer throughout this work to the two situations above as the the disjoint 
case and the non-disjoint case. Basically, the results in the disjoint case can be obtained by 
following the arguments in [1] while in the non-disjoint case the context is different because 
of the appearance of the Malliavin derivatives in the expression of square mean of (9). 

Proposition 2 Suppose H > | . Then for every r > the stochastic Skorohod integral in 
(16) is well-defined. 
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Proof: We will use the representation of (9) of Z^(t) as a Skorohod integral with respect 
to the measure W. Next we use the bound (11) with p = 2. Let us apply it to the process 

(s,r') ->• l [0) t](s)lA a (r) / dad 1 K(a,s)Q r (a) 

J s 

(which is a two-parameter process) and to H = L 2 (M x M + ,m). The variables t,r and the 
set ^2 with 7712(^2) = 1 are fixed. Below D is the Malliavin derivative with respect to the 
isonormal process W (see Section 2). We will get, by (11), 



E(Z*(t)) a < E f ds [ f dad 1 K(a,s)Q r (a^ 

JO Ja 2 V ) \Js 

da J —p- yD aj p J dadiK(a,s)Q r (a) 



'A 2 

rt r ft 
+E / ds 



= E^dsJ^ ^2 dadiK(a,s)Q r (a)^ 

f l f cdr' f* f cd/3 / f l 
+E / ds / — ^ / da / / dadiif(a, s)D a)l3 Q r {a 

Jo Ja 2 V ) Jo Ja 2 P \Js 
= E j Q ds j J^Jy da9iK(a,s)Q r (a)^ 
/"* f cdr' f l f cd(3 ( f l 

+ E / dfi / 77\2 \ da / ~tf2~ / da ^l- fi: ( a ' S ) ( 5r-(a)lc r (a)(a,/?) 

JO JA 2 l r ; Jo Ja 2 P \Js 
■= Ti+T 2 . 

since, by the chain rule of the Malliavin operator (13) 

D a> pQ r (a) = Q r (a)l Cria) (a, (3). (19) 
We need again to consider two cases. 

The disjoint case: We have A 2 f\ (0, 1) = (this happens when A 2 = (c, 00), c > 1 
for example). In this case the term denoted by T 2 vanishes because 

l Cr(a) (a,/3)U 2 (/3) = 0. 

We need to show that 



T\ = m 2 (A 2 )E J ds (j dad 1 K(a,s)Q 



2 

n a) 



is finite under condition (20) and this is exactly the computation in [1], proof of Proposition 
2.1. 
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The non-disjoint case: Suppose A 2 n (0, 1)^0 (this happens when A 2 = (c, 00) 
with c < 1). In this case we need also to show that 

T 2 < 00 

By bounding the indicator function l(j r ( a \(a, ft) by 1 (note that d\K(a,s) > for every 
a,s) 

E ( I dadiK(a, s)Q r (a)l Cr ( a )(u, P) 



< E (J dad 1 K(a,s)Q r (a] 



2 



and by computing the integrals dr' and d(3 

ft rt 



T 2 < m 2 (A 2 ) 2 E [ ds [ da ( ( dadiK (a, s)Q r (d 

JO JO \Js 

= m 2 {A 2 ) 2 E ds (J^ dad 1 K{a,s)Q r (a) 



= m 

and this is finite under (20). ■ 

4 The Multifractal Random Walk 

The purpose of this section is to study the limit as r — > of the family of stochastic processes 
(Z,?(t)) with fixed t > 0. This limit will be called the the Multifractal Random Walk. 

We will assume throughout this paragraph that Z^(t) is defined by (16) with 
given by (15). Moreover we will suppose 

m 2 (A 2 ) = 1 

and A 2 n (0, 1) / 0. The disjoint case A 2 n (0, 1) = follows from [1]. 
We will consider the following assumption: 

cip(2) + 2H > 1. (20) 

Remark 3 Since <p(2) = —1, the condition (20) means that c < 2H — 1. Since H > \ we 
can choose c > in order to have (20). By assuming (20), we are in the case (A) in [1]. 

We have the following limit theorem. 

Theorem 1 Assume (20). For every t > the sequence of stochastic process (Zf?(t)) r>0 
defined by (16) converges in L 2 (J7) to a random variable Z H (t). 
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Proof: Let us fix r, r' £ (0, 1) with r < r' . We will first compute the L 2 (fi) norm of the 
increment Z^(t) — Z% (t) where t > is fixed. We can write, with given by (15) 



E\Z?(t)-Z$(t)\ 



E 



E 



E 



(Q r (u)-Q r ,(u))dB* 

(Q r (a) - Q r ,{a))diK{a,u) ) dWj^ 



t / r t 



\Ju 



JA 2 \Ju 



< E du 



JA 2 



cdx 



da(Q r (a) — Q r '(a)) d\K{a,u) ) dW(u,x) 

2 



da (Qr(a) — Q r '(a)) d\K(a, u) 



E / du 



o 



cdx 



dv 



A 2 x Jo JA 2 V 



cdy 

/2 



D v , y / da (Q r (a) - Q r '{a)) d\K{a, u) 

J u 



where we used the bound (11) for the 1? norm of the divergence operator. Using the 
differentiation rule (19), we get 

E\Z*{t)- Z^{t)\ 2 



< m 2 {A 2 ) f du f da f dbd 1 K(a,u)d 1 K(b,u) 

JO J u J u 



xE(Q r (a) - Q r ,{a)) {Q r (b) - Q r ,(b)) 
rt r i r cdy r rt 



+m 2 (A 2 ) [ du f dv f 

JO Jo J A 
A(r,r') + B(r,r'). 



a 2 r 



dadiK(a,u) (Q r (a)l Cr (a){v, y) ~ Qr'{a)lc r ,( a ){v,y)) 



Let us first treat the term denoted by A(r,r'). Using property (10), 



t rt 



A(r,r') = m 2 (A 2 ) du da dbd 1 K(a,u)d 1 K(b,u) 

JO J u J u 

xE [Q r (a)Q r (b) - 2Q r (a)Q r {b) + Q r ,{a)Q r ,(b)] 
= m 2 (A 2 ) f du I da f dbdiK(a,u)d x K (b,u)E [Q r * (a)Q r > (b) - Q r (a)Q r (b)] 

Jo J u J u 

= eJ du^y dad\K(a, u)Q r i(a)^j du — E J du (^j dadiK(a,u)Q r (a)^j du 

= E(J*Qr*(u)dBA -E(J*Q r (u)dBA 

where B H denotes a fractional Brownian motion independent by W. The convergence of 
this term is similar to the study in [1] we are in the case (A) in their paper (see Remark 3). 
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The summand B(r,r') will be handled as follows. First, note that 

B(r,r') = m 2 (A 2 ) [ du [ dv [ ^ f da I dbd 1 K(a,u)d 1 K(b,u) 
Jo Jo Ja 2 y Ju Ju 



xE(Q r (a)l Cr(a) (v,y) - Q r >(a)l Cr , (a) (v, y)) (Qr(b)l Cr[b) {v,y) - Q r ,(b)l Crl{b) (v, 

m 2 (A 2 ) [ du f dv f ^ [ da f dbd 1 K(a,u)d 1 K(b,u) 
Jo Jo Ja 2 y Ju Ju 



IA 2 

xE [Qr{a)Qr(b)lc r (a)(.v,y)l Cr (b)(v,y) - Qr{a)Qr(b)lc r (a)(v,y)lc r ,(b)(v,y) 
-Qr(a)Q r (b)lc r ,(a){v,y)lc r (b){v,y) + Qr>(a)Q r >(b)lc r ,( a )(v,y)lc r ,(b){v,y)} 

where we used again (10). By decomposing 

C r (a) = C r ,{a) U (C r (a) \ C r ,(a)) 

we will have 

B(r,r') = m 2 {A 2 ) [ du [ dv [ ^ f da [ dbd 1 K(a,u)d 1 K{b,u) 
Jo Jo Ja 2 y Ju Ju 

xE [Q r , (a)Q r , (b) - Q r (a)Q r {b)\ l Cr , (a) (v, y)l Cr , (6) (v, y) 

+m 2 (A 2 ) [ du f dv [ ^ / da / dbd 1 K(a,u)diK{b,u) 
Jo Jo Ja 2 y Ju Ju 

><^Qr(a)Qr(b)lCr(a)\C r ,(a)(v,y)lc r (b)\C r ,(b)(v,y) 

:= Bi(r,r') + B 2 (r,r'). 

Obviously, by bounding the indicator functions by 1, 

ft rt r „j„, ft ft 



£i(r,r') < f du f dv f ^ f da [ dbdiK{a, u)d x K(b, 
Jo Jo J a 2 y Ju Ju 



>A 2 

xE [Q r ,{a)Q r >{b) - Qr{a)Q r {b)] 

and therefore it converges to zero as r, r' — > by using exactly the same argument as in the 
case of the term A(r, r'). 

Concerning B 2 (r, r'), since for < r < r' < 1 and y 6 A 2 

lc T (a)\C r ,(a)(v,y) < \r,r<)xA 2 (v,y) 
we can bound it in the following way 

B 2 (r,r') < m 2 {A 2 ) f du f dv j — [ da [ dbd\K(a, u)d\K(b, u) 
Jo Jo J a 2 y Ju Ju 

xEQ r (a)Q r (b)l {r y )xA2 (v,y) 
= m 2 (A 2 f(r' -r)E I du I da f dbd 1 K(a,u)d 1 K{b,u)EQ r (a)Q r (b) 

JO J u J u 
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and by interchanging the order of integration and using 




dud\K(a, u)d\K(b, u) = cn\a — b\ 



2H-2 



(21) 



and since 



EQ r (a)Q r (b) 



e -ip(2)m r (\a-b\) < e c^(2) log \a-b\ 



= \a — 



5|CV(2) 



we get 




dadb\a - b\ c ^ +: 



Jo Jo 

we clearly obtain than B2(r,r') goes to zero as r,r' 



under condition (20). 



Definition 1 The process (Z H ) t >o from Theorem 1 will be called as fractional Multifractal 
Random Walk. 



We will discuss here some immediate properties of the fractional MRW from Definition 1. 
Basically, as we mentioned before, in the disjoint case the fractional MRW Z H has the same 
properties as in the situation when Q and B H are independent: self-similarity, stationarity 
of increments and long-range dependence (see [1]). But in the non-disjoint case, we will 
show that even the 1? moment of the fractional MRW does not scale exactly. We provide 
an exact calculation in order to show this phenomenon. On the other hand, we can control 
the LP norm of the increment of the process and we find some kind of asymptotic scaling 
even in the non-disjoint case. 

5.1 Scaling of the second moment 

We first analyze the second moment of the increments of Z H . 
The disjoint case: Following [1], we have 



5 Properties of the MRW 
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where the meaning of the quantity EQo(a)Qo(b) is given by (5) with r = and where we 
used the identity (21). Then, for every h > 



pht pht 

E (Z H (ht)) = c H / db\a - b\ 2H - 2 E(Q (a)Q (b)) 
Jo Jo 

= h 2H c H [ [ db\a-b\ 2H - 2 E(Q Q (ha)Q (hb)) 
Jo Jo 

= h 2H Ee 2nh c H [ [ db\a - b\ 2H - 2 E(Q (a)Q (b)) 
Jo Jo 

= h 2H Ee 2Uh E {Z H {t)f = h 2H+c ^E {Z H (t)) 2 . 

where we used the scaling property (7). Actually, it is not difficult to see that for every 
p > 1 

E (Z H (ht)) P = h 2Hq Ee pQt = h 2H P+^P\ 

Moreover, we have the self-similarity (Z H (at)) t£ {Qi] a H+ ^ a (Z H (i))t e [o.i], the station- 
arity of the increments and the long-range dependence in the sense that EX^Xq ~ r 2 Hk 2H 
where X k = Z H ((k + l)r) - Z H (kr) with k > integer and r > 0. 

The non-disjoint case: 

The situation is different in the non-disjoint case and we will see that even the 
second moment of the fractional MRW does not scale. We can compute exactly the L 2 
norm of Z H (t). 



E\Z H (t)\ 2 = j ds (J daQ (a)diK(a,s)da 



x (v^P f s dbQ (b)d 1 K(b,s)db\ ^D 8<r J daQ (a)d 1 K(a,a)da 
f ds daQo(a)diK(a, s)da 



+ l >lsl ^l^Lw 



2 



x (^J daQ (a)d 1 K(a,a)l Co{a) {s,r) S j (J^ dbQ (b)l Co{b) {a, P)diK(b, s)db 
= I(t) + J(t) 

Consider h > 0. Then, applying the above formula to t = ht and making several changes 
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of variables, we will get 
E\Z H (ht)\ 2 = h 2H+c ^I(t) 

Jo Ja 2 r 2 Jo J A2 p 1 
x (^j daQ (a)diK(a,a)l Co (ha)(hs,r)\ (J^ dbQ (b)l Co ( hb) (ha, P)d x K{b, s)db ] 

The integral with respect to dp and dr can be computed explicitly. For example, when 
A2 = (c, 00), c < 1, we find 



cd/3. 



02 l c {hb)(ha,p) = c 
A 2 P 



1 



c V 2h\b - a\ 

and it is clear that the term J(t) does not scale exactly as I(t). 
5.2 The control of the increments 

Let us estimate the LP norm of the increment Z H (t) — Z (s) of the fractional Multifractal 
Random Walk introduced in Definition 1. Fix s,t £ [0, 1] with t > s. We will not insist on 
the disjoint case because the calculations in [J ] still hold, so we will have 

E \Z H (t) - Z H {s)\ P ~ C p \t - s \ 2Hp+c ^ 

with (p given by (6). 

Let us consider the non-disjoint case. In this case Z H (t), which can formally be 
written as J * Qo(y)dB H (u) is an anticipating (Skorohod integral). We need to use Meyer's 
inequalities (11) in order to estimate its LP norm. 

Actually 

Z H (t) - Z H (s) = f dWW(u)F s , t (u) 
Jo 

where we denoted by 

F 8tt (u) = l(o,t)(") / d 1 K(a,u)Q (a)da- l(o, a )(ti) / diK(a,u)Q Q (a)da. 

By Meyer's inequality (11) 

E\Z H (t)-Z H (s)\ P < E^duf ^-\F 8it (u)\ p 

Jo Ja 2 r 2 
1 r „a„ rl 



cdr f j f cd/3 
A(t,s) + B(t,s) 







2 



+ E / du I L da I. -f\Da,0F s , t (u) 
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The terms A is exactly the LP norm in the disjoint case, so 



A(t,s) < c p \t - s\ 2Hp+c ^ p) . 



Concerning the summand denoted by B{t, s), we proceed as in the proof of Proposition 2: 
we use first the Malliavin differentiation D a ^Q r (a) = Q r (a)lc (a, /3), then we bound the 
indicator function lc ( a > P) by 1, then we integrate dr and d/3 and we obtain 



because the right hand side is equal, modulo a constant, to A(t, s). Taking into account the 
above estimates, we conclude that 



for every t,s. 

6 Financial Statistics 

As mentioned in the Introduction, the multifractal random walks appears nowadays as a 
serious candidate to model the financial time series. In order to compare its behavior with 
real data, one needs to simulate it. The main difficulty consists in the fact that, in our 
construction, the variables Q and B H are dependent. From the theoretical point of view, 
the Malliavin calculus offers convenient techniques but thinks are complicated in from the 
practical point of view. 

6.1 Simulation Scheme 

Recall that the fractional MRW is defined by 



where the measure M and the fBm B are dependent. As throughout our pa- 
per, we choose M to be the Gaussian measure. Consider a random variable w such that 

e M (C(t)) _ (t) 
Ee M(C(t)) — e 

The first step consists in generating a sequence of random variables W = (w(ti))i 
such that w(ti-{-i) and w(U) are correlated and with Gaussian distribution A/" n (^, S) where 
the covariance matrix E depending on the auto-correlation function y(\t — (t + h)\) = 
Cov(w(t) , w(t + h)) . We consider the choice of 7 from [5] with j(k) decreasing with respect 
to k. Then, assuming that the matrix E is positively definite, there exists a triangular 
matrix C such that C*C = E. In this way, the components of 




E \Z H (t) - Z H (s)\ P < c p \t - s \ 2Hp+c ^ p) 



(22) 




(23) 
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W = fl+CWn(0,I), 



are correlated and with normal law. 

Let us explain the idea to simulate the integral (23). This integral is a divergence 
integral and in principle its simulation is difficult. But in our particular case, we can use 
the following approach. Let {0 = to < t\ < ■ ■ ■ < t n = t} a partition of the interval [0, t]. 
in the dependent case (meaning when W and B H are independent), (23) can be naturally 
approximated by 

n-1 

Y^e w ^{B H {t t+l )-B H (t t )). 

1=0 

Since the simulation of B H is well-known (see e.g. [6] for an explicit algorithm), we can 
generate the above sum. In the dependent case, using the integration by parts formula 
5(Fu) = F5(u) + (DF,u), 5 being the Skorohod integral with respect to B H (see [24]), the 
same sum can be expressed as 

n-1 

£ (e^\B H (t i+1 ) - B H { tl )) + (De«<M t l (ti ,t i+l) >) 
i=0 

n-1 

= J2 (e w{U) (B H (t l+1 ) - B H (t t )) + e w ^(Dw( tl ), l (u , u+l) )) 
i=0 

where the scalar product is in the Hilbert space associated to the fBm and, since w(tj) 
is Gaussian, Dw{ti) is deterministic and it depends on the set A<i (more exactly on the 
intersection of A2 and the interval (0,1) which is in principle very small). Thus, a possibility 
to simulate (23) is to approximate it by YaZq {e w{h) {B H {t i+ {) - B H {ti)) + ai e w ^) with 
suitable coefficients a;. 



6.2 Stylized Facts 

Let S(t) be the price at time t > T of an financial asset, X(t) = In S(t) the log price and 
then, log-returns at lag r are given by 



6 T X(t) = X(t) - X(t - r) = In 



S(t) 



S(t-n 

and, for sake of simplicity and with no loss of generality, we assume that log-returns are 
centered, E[5 T X(i)] = 0, thus, we can interpret the volatility as the squared log-returns. 

The construction of the proposed MRW allows us to verify several properties of 
stock market fluctuations, what we call stylized facts. The first two that we present are 
specific to the MRW, long memory in absolute log-return and in volatility; the others are 
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achieved thanks to our construction, leverage effect and long memory in return. 

The random walk that we introduce in our work allows the conditional variance and 
the random noise, which is a fBm with Hurst parameter H > |, to be dependent. This 
dependence should involve the so-called leverage effect, that is, the correlation between the 
log-return at time t and the squared log-return (the volatility) in the future. The case 
H > 1/2 implies that the noise B H exhibits long-range dependence. This should involve a 
long-range dependence of the fluctuations of the financial time series. These facts (leverage 
effect and long-range dependence) will be checked on the data. Actually the leverage effect 
is empirically observed on financial index such as the SSzP 500 (US market), Nikkei 225 
(Japan market), FTSE 100 (UK market), CAC 40 (french market) etc. But it does not 
appear on the Forex market and on the commodities, for example. We also notice the 
absence of auto-correlation in log-returns if we consider low frequency financial data but 
this auto-correlation structure clearly appears in the case of square or absolute log-returns. 
Nevertheless, at high frequency, we observe a dependence relation which is going to be 
formed even for usual (not squared and not absolute) auto-correlation function, (see e.g. 
[19] or [25]). 

The data that we used in order to check the leverage effect and the long-range de- 
pendence come from S&P 500 index with a frequency of 15 seconds from 2012-02-28 to 
2012-06-26, 131011 points. We compare the characteristics of the S&P 500 index with the 
simulation of our multifractal random walk. As we mention, we compare our simulations 
with the data from the S&P 500 index. In fact, we only present the comparison for the 
auto-correlation function since the multifractal structure of the asset return distribution 
are already well known and well describe by standard multifractal random walk (see e.g. 
[22]). The numerical results are presented in Figures 1 and 2. To compare, we also present 
the simulation of the fractional MRW with H = 0.62. This value has been obtained by a 
classical R/S test. 

Since the Hurst parameter is bigger that 0.5, we should have a long-range dependence 
for the log-returns. The auto-correlation function of the log-returns is defined by 

C(k) = Coir(5 T X(t),5 T X(t + k)). 
Usually, we say that the log- return have long-range dependence if it satisfies, 

C(k) ~ k~ a , a > 0, 

We give in Figure 3 the empirical auto-correlation function of the log-returns, the best fit 
obtained in power law and the simulated auto-correlation function of the fractional MRW 
with H = 0.62. 

We also notice a long-range dependence for the log-returns in the case of S&P 500 in- 
dex at high frequency, but this less than for the squared or absolute long-returns. We could 
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V 

20000 40000 60000 80000 1 00000 1 20000 

time (15 seconds) 

Figure 1: S&P 500 index , at 15 seconds, from 2012-02-28 to 2012-06-26, 131011 points. 

advance two main ideas to explain this small persistence in high frequency log-returns. 
First one is that all agents include information, good or bad news, in their investment 
strategies, but there may have some latencies in the response of investor and cause this 
auto-correlation. The second one can be cause by huge investors such as mutual fund who 
want to invest significant sums and place large order. They have to split their large order 
in a sequence of smaller orders to acquire the number of shares wanted, and then, could 
cause some persistence. A simple and current strategy to do that is called iceberg order, 
more sophisticate could be found in the theory of 'optimal trading'. 

The persistence of dependence for the log-returns it not easily observed, but the 
presence of the cluster of log-returns indicates basically that strong variations are followed 
by other strong variations. Then, we should observe a persistence in absolute and squares 
log-returns. As before, we present in Figure 5 and Figure 4 the empirical results for SSzP 
500 with its fit in power law and the simulation of the MRW {H still equal to 0.62). 

We now present the second result that we obtain via Malliavin calculus construction. 
We have assumed that the measure M and the integrator B H was dependent, hence, the 
MRW should verify the correlation between return at time t and squared return (volatility) 
around t, say t + k<t'<t + k, this quantity correspond to the leverage effect, 

C(t,t') =E [5 T X(t)(5 T X(t')) 2 } , t + k<t'<t + k, (24) 

This quantity is usually normalize by E[(5 r X(t)) 2 ] 2 , we use this normalized formulation, 
see [11] for details. In the literature, one can usually find studies of the stylized facts at 
low frequency, based on financial indices or on the average of financial indices (see [10] for 
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Figure 2: Top: log-return of S&P 500, at 15 seconds, from 2012-02-28 to 2012-06-26, 
131011 points. Down: simulation of the MRW with H = 0.62. 

a very interesting empirical study), that's why we have chosen to work with the S&P 500. 
We will follow the same approach, but we also present results for high frequency. 

The empirical study of the leverage effect is given in Figure 6. To compare, we also 
plot the leverage effect for the fractional MRW in Figure 6 at low frequency. 
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k (15 seconds) 



Figure 3: Auto-correlation of log-returns, k from 1 to 100. In black the S&P 500, in grey 
the confidence interval, in red the MRW and in blue to fit in power law (a = 2.1). Notice 
that the mean square error between the data and the numerical results is 3.868 x 10 -4 and 
if we take into account the fit, it is 0.404 x 10 



III! Illilliilil 

llifjliiiHIlll] 



v° 



TtTTtTttTttT 



>v».*» 



WW 



>V,\°/c° "A-.s 

WW 



k (15 seconds) 



Figure 4: Auto-correlation of the absolute log-returns, k from 1 to 100. In black S&P 500, 
in grey the confidence interval, in red the MRW and in blue the fit in power law (a = 1.5). 
The mean square error is 1.720 x 10~ 4 and with the fit it is 0.223 x 10 -4 
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20 40 60 80 100 

k (15 seconds) 



Figure 5: Auto-correlation of the squared log-returns, k from 1 to 100. In black S&P 500, 
in grey the confidence interval, in red the MRW and in blue the fit in power law (a = 2.4). 
The mean square error is 2.834 x 10 -4 and with the fit 1.209 x 10~ 2 
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